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THE NATURE OF THE GOULD BELT FROM A 
FRACTAL ANALYSIS OF ITS STELLAR POPULATION 
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ABSTRACT 

The Gould Belt (GB) is a system of gas and young, bright stars distributed 
along a plane that is inclined with respect to the main plane of the Milky Way. 
Observational evidence suggests that the GB is our closest star formation com- 
plex, but its true nature and origin remain rather controversial. In this work we 
analyze the fractal structure of the stellar component of the GB. In order to do 
this, we tailor and apply an algorithm that estimates the fractal dimension in a 
precise and accurate way, avoiding both boundary and small data set problems. 
We find that early OB stars (of spectral types earlier than B4) in the GB have 
a fractal dimension very similar to that of the gas clouds in our Galaxy. On the 
contrary, stars in the GB of later spectral types show a larger fractal dimension, 
similar to that found for OB stars of both age groups in the local Galactic disk 
(LGD). This result seems to indicate that while the younger OB stars in the GB 
preserve the memory of the spatial structure of the cloud where they were born, 
older stars are distributed following a similar morphology as that found for the 
LGD stars. The possible causes for these differences are discussed. 
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INTRODUCTION 



The brightest stars near the Sun are mainly distributed along two great circles over 
the sky: the Milky Way plane and another strip of bright stars known as the Gould Bel t 
(GB) which is inclined respect to the Galactic plane (lGouldlll879l : IStothers &: Frogellll974j ). 
Alt hough th e exist ence of the GB was formerly reported as soon as in 1847 by Sir John F. 
W. iHerschell (118471 ) from naked eye observations of the southern sky, its nature and origin 
are still poorly understood. The GB is a complex system of gas and stars composed not 
only by single massive stars but also by OB associations and an interstellar medium (dust, 
neutral hydrogen, and molecular clouds) that shows some kinetic fea tures which have bee n 
considered by several authors as connected to the stellar population (ILindblad et al.lll973l ). 
The GB was firstly interpreted as an apparently expanding ring (or torus) of young massive 
stars, but later, with the detection of yo ung low-mass stars, i ts structure seemed to be better 
represented by an inhomogeneous disk (iGuillout et al.lll998l ). Recent studies characterize it 
as a stellar disk having elliptical shape, with a major semiaxis of ~ 6 00 pc and rainor se miaxis 
of ~ 400 pc, which is inclined around 18° to the Galactic plane (lElias et al.ll2006a). Th e 



stellar component exhibits a range of ages lesser than 100 million years (iTorra et al.ll2000l ). 
Given its size (~ 1 kpc), age range, and relationship with the gas in differ ent phases, it has 



been suggested that the GB would be our closest s t ar for mation complex (lElmegreen et al. 



20001 ). in the sense given to this concept by lEfremovl (119951 ). i.e., the largest region of a galaxy 
showing a coherent star formation process, where the term "coherent" should be interpreted 
as originated from the same, monoparental, giant gas cloud. 

However, there exist both theoretical and observational results showing inconsistencies 
with this scenario: the age pattern of OB a ssociations is irregular, with the youngest ones 
located well outside the expanding gas ring (jPerrot fc Grenierl |2003| ) : the kinematics of the 
gas does not seem to be in agreement with the velocity field of the stellar component; 
and the modelization of gas dynamics is not compa tible with the estimated radius of the 



stella r component, as well as with its range of ages (jPerrot fc Grenien l2003l : iMoreno et al. 



19991 ). Moreover, the analysis of t he stellar spatial velocity diagram in the U ~ V plane 
shows a clear bimodal distribution (lElias et al.ll2006bl ). All these results make it difficult to 
imagine a single monoparental origin for the GB and suggest a GB formed by a spurious 
concatenation of different stellar subsystems whose formation processes were not necessarily 
connected among them. This point is particularly important for a full understanding of 
the mechanisms involved in the formation of the Milky Way. Can we design some kind of 
experiment able to answer this question? Here we propose to analyze the fractal structure of 
the stellar component of the GB in order to go deeper into its nature and origin. The base of 
this approach lies in the fact that the interstellar medium see ms to show a fractal structure 
when observed at the scale of clouds (IBergin fc Tafallal 120071 . and references therein) and a 
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multi-scaling behavior (jChappell &: Scald I2OOII ) when galactic scales are considered. Thus, 
stars forming from the same cloud s hould exhibit fractal patterns t oo if their birth places 
uniformly follow the densest regions (lElmegreen &: Elmegreen Il200ll ). In this work, we first 
develop a "reliable" algorithm to compute fractal dimensions of a sample of discrete points 
(§ [2]) and then we use it to study the distribution of stars in the GB (§ [3]). The main 
conclusions are summarized in § |H 



2. ESTIMATING THE FRACTAL DIMENSION 

Strictly speaking, a fractal is defined as an ob ject whose Hausdorff dimension is larger 



than its topological dimension (IMandelbrotl 1 19831 ) . To estimate the Hausdorff dimension 



authors use different working definitions that fit their methods and needs, and thus there is 
not a unique definition of fractal dimension. When dealing with a distribution of points i n 



space it is very useful the so-called "correlation dimension" (IGrassberger fc Procaccialll983l ). 
This is a widely used method because of its robustness and because it is relatively easy to 
implement on real (experimental or observational) data. 

Let us consider a distribution of points in space with positions x. The number of 
other points within a sphere of radius r centered on the i-th point is given by the expression: 



N 

j=l, j^i 



where H{x) is the Heaviside step function. This number can be evaluated by choosing M 
different points as centers and then averaging to obtain the probability of finding a point 
within a sphere of radius r centered on another point. This probability is expressed in the 
form 

1 ^ 

For a fractal set this quantity, called the correlation integral, scales at small r as 

C(r) ~ , (3) 

being Dc the correlation dimension, which in practice can be identified with the slope of the 
best fit in a log C(r) — logr plot. In other words, the correlation dimension tells us how 
it changes (as r increases) the probability that two points chosen randomly are separated 
by a distance smaller than r. For a homogeneous distribution of points in space we expect 
Dc = 3, whereas in a plane Dc = 2. If points are distributed obeying a fractal geometry then 
Dc < 3 (in 3D-space) or Dc < 2 (in a plane). 
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Typically, when evaluating Dc for real data (and not for infinite, perfect fracta ls) the 
powe r-law scaling relation (eq. [3]) is valid only within a limited range of r values ((Smith 



19881 ). At very small scales (of the order of the mean distance to the nearest neighbor) the 



distribution looks like a set of isolated points and D^. tends to zero. On the other hand, the 
finite size of the set also results in decreased values at large r values (of the order of the 
set size), because near the edge each point is surrounded by other points only on one side and 
then C{r) tends to be underestimated. It has been proved that these and many other effects 



can lead to bad estimations of Dr., mainly w hen too few data are available (jSmithI Il988 



Kitoh et al.l 120001 : ICiccotti &: Mulargial l2002l ). In order to assess how the limited number 



of data points would alter our estimation of Dc for the GB, we have done some numerical 
experiments. Figure [1] shows log(C) — log(r) plots for random distributions of points within 
disk-like structures. Boundary effects were avoided by keeping the condition that sampling 
spheres always are inside the volume occupied by the disk. For the three-dimensional case 
(Fig. [1^), the expected linear behavior is observed at high r values for a relatively high 
number of points A^, but at low r values we s ee departure s from the linear behavior arising 
from the lack of statistic of finite samples JSmithlEoSsI ). This tendence becomes more 
evident as decreases because the mean distance between nearest neighbor increases. For 
the case A^ = 200 points it is very difficult to infer a linear behavior at all. This is an 
important problem because many times the number of available data points is rather small. 
In these cases the estimated fractal dimension de pends strongly on the spatial range used 
to calculate the slope in the log(C) — log(r) plot (jPisarenko fc Pisarenkd 119971 ). However, 
when the calculation is done on the projected distributions (Fig. [T}d), the range of r values 
used to estimate the integral correlation can be extended to higher values, while still keeping 
the spheres (that actually become circles in the 2D projection) inside the sample. In this 
case, the linear behavior is clearly appreciated -at relatively high r values- even when only 
A^ = 200 points are used. The average D^. values for ten random realizations, calculated by 
doing linear fits in the range logr > 1.5, are shown in Tabled! The results are always close 
to the theoretical values (3 or 2) but the standard deviations become quite high in the 3D 
case for a relatively low number of points. In any case, the "problematic" region of relatively 
low r values has to be excluded of the linear fit. In order to do this, we impose the condition 
that the standard deviation of each C(r) value must be smaller than the corresponding C(r) 
value. This simple criterion eliminates poorly estimated data (i.e., bad sampling) ocurring 
mainly at too low r values. 

As an example. Figure [2] shows a random distribution of 100 points in a square region of 
side 1000. In order to take into account edge effects, we first determine which points of the 
sample are vertices of the convex hull, i.e. we determine the minimum-area convex polygon 
containing the whole set of data points (square symbols connected by lines in Figure [2]). For 
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this we use the algorithm proposed by lEddyl (Il977l ). Now we can place circles of different 



radii r to evaluate C{r) according to equation ([2]). The result is plotted in Figure [31 The cal- 
culations were carried out both taking (open circles) and not taking (asterisks) into account 
boundary effects. If circles are allowed to cross boundaries C (r) tends to be underestimated, 
and this effects is higher the higher the radius r. For this reason Dc would tend to be smaller 
than the expected value if this effect is not considered (for comparison, solid line in this figure 
shows the expected = 2 result). The range of r values that fulfill the condition that the 
standard deviation of C{r) is smaller than C{r) is indicated in Figure [3] as vertical arrows. 
For r values below this range the departures from the linear behavior become more evident. 

The algorithm we have developed follows the next steps: first, the three-dimensional 
distribution of points is projected on its mean plane; the boundary is determined finding 
the convex hull of the sample; the correlation C(r) is calculated using equations ([I]) and 
([2]) but always keeping the circles inside the sample boundary; and then the correlation 
dimension is determined as the slope of the best logC(r) — logr linear fit. Very low r values 
(and consequently poorly estimated C(r) values) are excluded from the fit. Finally, we use 
bootstrap techniques to estimate the uncertainty of the calculated value: we repeat the 
calculation of on a series of random resamplings of the data and the standard deviation 
of the obtained set of fractal dimensions is taken as the error in our estimation. When 
working with relatively thin disks the situation is practically the same as taking a two- 
dimensional slice of a spherical, three-dimensional distribution of points within a sphere. In 
this case, the fractal dimension of the three-dimensional distribution of points D JSD) and 



the t wo-dimensional dimension Z)c(2D) are related through the simple expresion ([Falconer 



1990[): 

D,(2D) = D,(3D) - 1 . (4) 

We have simulated two-dimensional fractal distributions of points in order to test the algo- 
rithm performance. The fractals were simulated by placing four squares of radius R/L (with 
L > 2) inside a square of radius R (we place one square in each quadrant). The procedure 
is repeated successively 10 times to obtain 4^° (~ 10^) points distributed according to a 
fractal pattern with dimension Dc = log 4/ logL. Finally, we randomly removed points from 
the fractal until reaching a given sample size A^. Figure [1[ shows some example results for 
the cases = 1.5 (L ^ 2.5) and D^. = 2 [L = 2). Even for relatively small sample sizes 
(A^ ~ 100) the algorithm is able to estimate Dc in a reliable way. 
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THE FRACTAL DIMENSION OF THE GOULD BELT 



The first problem th at arises wh e n stud ying; the Gould Belt is how to select the stars 
belonging to this system. lElias et al.l (l2006al ) developed a new method to perform a three- 
dimensional spatial classification that was applied to a sample of 553 OB stars from the 
Hipparcos catalog with precise distances of less than 1 kpc. This allowed them to separate 
and estimate the spatial structure of the stars belonging to t he GB and to the LGD. The 
distributions of MK spectral types for the sample of stars used ( jElias et al.ll2006al ) are shown 
in Figure [51 We see two populations clearly differentiated, one with spectral types earlier 
than B4 (which we will call "early" stars) and another one with spectral types B4-B6 ( "late" 
stars), that can be associated with two age groups centered at ~ 20 Myr and ~ 70 Myr, 
respectively. The spatial distribution of these stars is shown in Figure [61 The projections on 
each of the mean planes clearly reveal a clumpy, filamentary structure for the GB (Fig. [6^), 
whereas the stars in the LGD seem to be distributed more homogeneously (Fig. [6]d). The 
almost free of stars gap seen in Figure M corresponds to the line of nodes in which both the 
GB and the LGD coexist. Close to this line, the probability of belonging to the GB has 
a high uncertainty because the bayesian probability in this region depends only on the "a 
priori" probability (i.e., the relative fre quencies of both sy stems) and not on the probability 
conditioned to the spatial position (see lElias et al.ll2006al ). This gap is almost imperceptible 
for the GB and notorious for the the LGD simply because in the sample of stars used there 
are more stars in the GB than in the LGD. What we have done in this work is to quan- 
tify the degree of inhomogeneity by calculating the fractal dimensioiil] for the GB and LGD 
(both early and late stars). The results are shown in Figure [71 For all cases the behavior is 
almost perfectly linear in this log-log plot, with correlation coefficients ~ 0.99. The number 
of circles used to evaluate C{r) for each r (i.e. M in eq. [2]) is usually of the order of 10^. 
However, strictly speaking, M depends on r in the sense that larger r allows fewer circles 
within the borders. For the largest r values considered M is of the order of 10. The slopes 
of the linear fits give the following fractal dimension values: 



GB-early: D^W) = 2.68 ± 0.04, 

GB-late: D^W) = 2.85 ± 0.04, 

LGD-early: D^{3D) = 2.89 ± 0.06, 

LGD-late: D^{3D) = 2.84 ± 0.06. 

If the empty gap mentioned before were producing some bias in the determination of 



^The term "fractal dimension" alone may be ambiguous because there exists several different definitions 
of this quantity. We will indifferently use the term "fractal dimension" or "correlation dimension" to refer 
to Dc- 
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for the LGD then the unbiased values should be even higher than the obtained ones. 



We see that the distribution of stars in the solar neighborhood exhibits a certain degree 
of fractality, with 2.7 < Df < 2.9. Statistical tests show that there is a difference between 
early stars in the GB and the other subsets that is, in the worst of the cases, statistically 
significant. Therefore, interestingly, the fractal dimension of early stars in the GB (-D/ ~ 2.7) 
is significantly smaller than that of the rest of the sample. Thus, these early stars have similar 
ages and are distributed following fra ctal patterns analogous to those observed in the gas 
of star forming regions in the Galaxy (jSanchez et al.ll2005l . 120071 ). Therefore, it seems very 
likely that this group of stars was born from the same parent cloud. In contrast, later stars 
in the GB have a somewhat different, more homogeneous type of distribution with a fractal 
dimension similar to that obtained for the LGD {Df > 2.85). 

What is the origin of this difference? First, we have to point out that, in principle, 
different global properties are expected for the interstellar medium (ISM) at different spa- 
tial and/or temporal scales because the main physical mechanisms acting at each scale 
are not necessarily the same. A monolithic gas cloud can be character ized as a turbulent 



medium, and this turbulence can be driven by many energy sources (lElmegreen fc Scalo 



2004 ). However, at larger spatial scales, the ISM in a spiral arm is influenced by other 
disturbances (gravitational, magnetic, etc.) that modify the gas distribution. The fractal 
analysis gives us a simple but objective measurement of such structure through the degree 
of smoothness or clumpyness. This kind of analysis has been applied extensively to evaluate 
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amples]. But often the analysis techniques are so diverse that it is not quite easy to extract 
robust co nclusions. It seems that th e average fractal dimension of single gas clouds is around 
Df ~ 2.7 (ISanchez et al.ll2005l . 120071 ). although it may vary when the spati al scale is increased 



and when the distribution of gas clouds in the galactic disk is analyzed (IChappell fc Scalo 



20011 ). But what about the stars? Young, new-born stars will reflect the conditions of the 
ISM from which they were formed. Therefore, a group of stars born from the same mono- 
lithic cloud, i.e. born at almost the same place and time, should have a fractal dimension 
similar to that of the parent cloud. This is exactly what we find for early-type stars in the 
GB: we obtain a fractal dimension very similar to the ISM value. Otherwise, if the star 
sample is representative of a population born from various clouds and/or with different star 
formation histories, then the fractal dimension will represent the gas distri bution at different 



spati al or temporal scales according to a multi-scale structure of the ISM (IChappell &: Scalo 



20011 ). At galactic level, the fractal dimension of the distribution of stars and/or star forming 
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regions exhibits a very wide range of values, but until now no correlation ha s been clearly 



found between fractal patterns at galactic level and other galactic properties (lOdekonI 12006 



Feitzinger fc Galinskil 119871 : iParodi &: Binggelill2003l ). There are, however, some suggestions 



time after the star formation process (Ide la Fuente Marcos & de la Fuente Marcos 


2006a 


b 


3 


Odekon 


2006; 


Schmeja & Klessen 


2006). probablv due to the action of some phvsical mech- 



anisms which tend to reorganize/destroy the original structure. 

Here we have found quite different values for the fractal dimension of early-type stars 
in the GB and the LCD. This difference reflects the fact that most of OB associations in the 



solar neighborhood are located in the GB (Ide Zeeuw et al.lll999l : lElias et al.ll2006al ). This 
is, the distribution of OB stars in the LGD does not seem to show any kind of "typical" 
stellar grouping, whereas young OB stars in the GB exhibit a more hierarchized distribution 
with multiple OB associations. In this sense, we can say that the GB shows clear signatures 
of the internal structure of its parental cloud. This clustering in the young population of 
the GB is quantified through the fractal dimension, whose value is very close to the value 
characterizing gas clouds in the Galaxy. Thus, the possible causes for different values for the 
fractal dimension of early-type stars in the GB and the LGD can be studied by analyzing the 
possible causes for different number of associations in both stellar systems. There are three 
possible factors (or a combination of them) that can contribute to this difference. Firstly, 
the age difference between both groups could be significant enough to explain this result. It 
is clear that the separation according to the spectral type of the stars provides only a gross 
age separation. The spectral types allow to derive an upper limit for the age but they do 
not give the real lifetime of a star that also depends on the initial chemical composition. 
Although the earlier spectral groups of both the GB and the LGD have similar mean values 
and variances (see Figure [5]), there exists the possibihty that one group is actually younger 
than the other one. This age difference could be enough to allow the LGD group to disperse 
showing a less structured distribution, whereas the GB group would still trace the density 
peaks of the parental cloud. Secondly, the difference in the degree of clustering among both 
groups can be produced by different physical conditions of the parental clouds. Much of 
stellar dispersal from active star-forming regions occurs on timescales of 10^ — 10^ years, and 
there are se veral possible scenarios to a ccount for this dispersal of young stars from their 
birth places (jMamajek fc Feigelsonll200ll ). The time scale for dispersal depends on the range 
of values of the densi ty peaks inside th e cloud which, at the same time, should depend on the 
local ISM pressure (jElmegreeru |2006| ) . Thus, clouds having different dynamical evolutions 
in different environmental conditions would likely give rise to stellar groups with different 
lifetimes. Depending on these critical values and on the age range of the stellar sample, 
we can get either a complex of associations or a uniform distribution of stars. Thirdly, 
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different gas distributions in the parental clouds can be responsible for variations in the 
fractal dimension among young stars in the GB and the LGD. For the GB, Dj agrees well 
with the range of v alues expected for the case of single monoparental clouds in the Galaxy, 
i.e., -Dj ~ 2.5 — 2.7 (ISanchez et al.ll2005l . 120071 ). On the other hand, the fact that the fractal 



dimension of the LGD is higher than the expected for individual clouds suggests that this 
population did not have a monoparental origin. Even though the young stars in the LGD 
were born in individual clouds each one having internal fractal dimensions around ~ 2.6, the 
set of clouds that formed this group were distributed following spatial pat tern with a higher 



fract al dimension at a larger scale according to a multi-scale scenario (IChappell &: Scalo 



20011 ). In other words, even though the two young stellar groups in the GB and the LGD 
were born at the same time, the structure of the parental gas at the birth time has likely 
driven the final geometric structure of these stellar populations. 

In order to restrict even more the problem under consideration, we have selected sub- 
samples of stars with spectral types earlier than B2 both for the GB and for the LGD. Stars 
belonging to these groups have ages younger than ~ 10 Myr. The determination of the frac- 
tal dimension for these "very early" subsamples yields similar results: Dc{3D) = 2.62 ± 0.08 
for the GB and Dc{3D) = 2.84 ±0.08 for the LGD. The uncertainties are higher than before 
because of the smaller number of stars in the samples. Obviously this result does not rule 
out the possibility that both groups have different ages but, in this case, it is more difficult 
to interpret these results in terms of age differences among groups. 

With the available data we cannot unambiguously discriminate the most likely expla- 
nation for the different fractal dimensions found for the very early-type stars in the GB and 
the LGD. Maybe all of these factors are contributing in some way to the problem. However, 
we consider that the main factors are probably related to different physical conditions in 
the parental clouds or different internal structure for the gas associated with a multi-scale 
scenario. Finally, we want to ment ion that besides t he different geometry and kinematics 
shown by the GB and the LGD (see EhaseLaDlsQoSfl, and references therein) the internal 



distribution of stars with spectral types earlier than B4 (quantified through the fractal di- 
mension) allows us clearly to differentiate both systems. This difference could be explained 
in terms of a hierarchical star formation scenario. 



4. CONCLUSIONS 

We have developed an algorithm which uses equations ([1]) and ([2]) to estimate the 
correlation dimension of the stellar component of the GB and the LGD. The novelty of 
the algorithm lies in the implementation of objective criteria to avoid boundary effects and 
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finite data problems at small scales. We find that early OB stars (of spectral types earlier 
than B4) in the GB have a fractal dimension ~ 2.7. This values is very similar to that of 
the gas clouds in our Galaxy. This result seems to indicate that younger OB stars in the 
GB "preserve the memory" of the spatial structure of the cloud where they were born. On 
the contrary, stars in the GB of later spectral types show a larger fractal dimension, very 
similar to the value found for stars of both age groups in the LGD (~ 2.8 — 2.9). Several 
factors (or a combination of them) can contribute to these morphological differences: age 
differences among the samples, different environmental conditions in the birth places, or 
different internal structure of the gas in the parental clouds. 

We want to thank the referee for his/her valuable comments and criticisms that improved 
this paper. E. J. A. acknowledges funding from MEG of Spain through grants AYA2004- 
05395 and AYA2004-08260-G03-02, and from Gonsejeria de Educacion y Giencia (Junta de 
Andalucia) through TIG- 101. 
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Table 1. The fractal dimension for random point distributions. The results are the 
average of ten random realizations. Here we show (left to right columns) the number of 
points in each test distribution {N), the correlation dimension (Dc) with its standard 
deviation (a) calculated both in the three-dimensional space (3D) and in the projected 

two-dimensional space (2D). 







a(3D) 


A:(2D) 


f7(2D) 


2000 


3.04 


0.10 


2.00 


0.01 


1000 


3.05 


0.27 


2.00 


0.03 


500 


3.13 


0.59 


1.98 


0.04 


200 


3.12 


1.32 


2.08 


0.14 
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Fig. 1. — Correlation integral C(r) for points distributed randomly within a disk of radius 
Rd = 2500 pc and half-heigh Zmax = 100 pc. The number of points is = 2000 (squares), 
1000 (circles), 500 (triangles), and 200 (rhombuses). All the curves have been arbitrarily 
shifted downward (except the upper one) for clarity, (a) The calculation is done in the three- 
dimensional space using spheres with radius ranging from the minimal distance between 
two points to a maximum of r^ax — -^max- As a reference, the dashed lines show the 
expected slopes of 3 for these cases. For the case A?" = 200 points, it is difficult to infer the 
expected linear behavior. Panel (b) shows the results when the calculation is done over the 
distributions projected on the Z = plane and using circles with different radii r up to a 
maximum of r^ax = Rd/'^- In these cases the expected value of the slopes (dashed lines) is 
2. The linear behaviors are clearly observed even for a relatively low number of points. 
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Fig. 2. — Example of random distribution of 100 points in a square region of side 1000. The 
vertices of the convex hull (see text) are indicated by squares connected by lines. The circles 
(with radius 200) illustrate how the sampling is done by keeping them inside the boundary 
defined by the convex hull. 




- 16 - 




1 1.5 2 2.5 3 

log r 

Fig. 3. — Correlation integral C{r) for the same point set shown in Figure El when the 
calculations are performed by taking (open circles) and not taking (asterisks) into account 
boundary effects. The solid line indicates the expected slope for a random distribution of 
points {Dc = 2), and the vertical arrows indicate the range for which the standard deviation 
of C(r) is smaller than C{r). 
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Fig. 4. — Correlation dimension as a function of the sample size N for fractal distributions 
of points. Each result is the average of 10 different realizations, and the bars are the average 
of the uncertinties calculated using bootstrapping (see text). The horizontal lines indicate 
the fractal dimension used to generate the distribution of points (1.5 and 2). 
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Fig. 5. — Distribution of spectral types for the sample of stars taken from Elias et al. f 2006 al). 
Continuous line refers to stars belonging to the GB, whereas dashed line refers to the LGD. 
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Fig. 6. — Positions of the stars in our sample. Stars are represented as circles (for the GB) or 
squares (for the LGD). The blue symbols refer to the spectral range 0-B3 and red symbols to 
the range B4-B6. The panels show the projections on the mean planes for each distribution, 
both for the GB (a) and for the LGD (b). 
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Fig. 7. — Projected correlation integral for stars in the solar neighborhood. The calculations 
were done for GB early stars (circles), GB late stars (squares), LGD early stars (triangles) 
and LGD late stars (rhombuses). Filled symbols indicate the range used for the linear fits 
(shown as lines). All the curves have been arbitrarily shifted downward (except the upper 
one) for clarity. 



